Magnetic Field Effects on the 
Superconducting and Quantum Critical 
Properties of Layered Systems with Dirac 

Electrons 

E. C. Marino and Lizardo H. C. M. Nunes 

Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
Cx. P. 68528, Rio de Janeiro- R J 21941-972, Brazil 

Abstract 

We study the effects of an external magnetic field on the superconducting prop- 
erties of a quasi-two-dimensional system of Dirac electrons at an arbitrary tempera- 
ture. An explicit expression for the superconducting gap is obtained as a function of 
temperature, magnetic field and coupling parameter (Ar). Prom this, we extract the 
B x Ar, T x Ar and B xT phase diagrams. The last one shows a linear decay of the 
critical field for small values thereof, which is similar to the behavior observed ex- 
perimentally in the copper doped dichalcogenide Cu x TiSe2 and also in intercalated 
graphite. The second one, presents a coupling dependent critical temperature T c 
that resembles the one observed in high-T c cuprates in the underdoped region and 
also in Cu x TiSe2- The first one, exhibits a quantum phase transition connecting a 
normal and a superconducting phase, occurring at a critical line that corresponds 
to a magnetic field dependent critical coupling parameter. This should be observed 
in planar materials containing Dirac electrons, such as Cu x TiSe2- 
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1 Introduction 



A lot of attention has been devoted recently to quasi-two-dimensional con- 
densed matter systems presenting a band structure such that the dispersion 
relation of the active electrons corresponds to the one of a relativistic massless 
particle. The kinematics of such electrons is described by a Dirac instead of a 
Schrodinger term in the hamiltonian [1,2]. This fact has a profound impact on 
the physical properties of the system. The Fermi surface reduces to a point, 
the Fermi point, where the density of states vanishes. This will drastically 
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affect the physical properties of such materials. In a previous work [3] , we in- 
vestigated the superconducting properties of a quasi-two-dimensional system 
of Dirac electrons and showed that they are completely different from the ones 
presented by usual Schrodinger electrons. There is, in particular, a quantum 
phase transition connecting the normal and superconducting phases, which is 
controlled by the magnitude of the effective superconducting interaction cou- 
pling parameter. The Cooper theorem, therefore, is no longer valid, as one 
should expect in the absence of a Fermi surface. 

Among the materials presenting Dirac electrons as their elementary excita- 
tions, there are a few, which have been intensely focused lately. These are 
high-Tc cuprates [4,5,6,7], graphene [8,9,10], carbon nanotubes [11] and tran- 
sition metal dichalcogenides [12,13]. 

The interplay between superconductivity and magnetism is a subject of central 
interest in any research involving superconducting materials. In particular, a 
key issue, both from the basic and applied physics points of view is the analysis 
of the effects of an external magnetic field on the superconducting properties 
of a system. In the present work, we study the effects of an applied constant 
magnetic field, perpendicular to a quasi-two-dimensional superconducting sys- 
tem containing Dirac electrons. As investigation method, we use the effective 
potential for the superconducting order parameter, which we evaluate both at 
T = and T ^ 0, as a function of the applied field (B). 

We firstly consider the zero temperature case and explicitly obtain the super- 
conducting gap as a function of the magnetic field and of the renormalized 
effective physical coupling parameter (Ar) that controls the superconducting 
interaction. This allows us to obtain the (T = 0) (B x Ar) phase diagram, 
which presents a quantum critical line separating the normal and supercon- 
ducting phases. A renormalization group analysis is then performed, demon- 
strating that the physical results do not depend on the renormalization point. 

We then study the nonzero temperature regime and explicitly obtain the su- 
perconducting gap as a function of the temperature and of the applied field. 
From this, we extract an implicit equation relating the critical temperature 
with the critical magnetic field, T C (T C , B c ). This allows us to obtain the T x A R 
phase diagram, which shows the critical temperature line as a function of the 
coupling, for different values of the magnetic field. This line resembles the one 
observed in the underdoped regime of high-T c cuprates and also in transition 
metal dichalcogenides [14] 

Also from T C (T C , B c ), we can derive the B x T phase diagram, for different val- 
ues of the coupling parameter Ar. The critical line separating the normal and 
superconducting phases, in this case, shows a linear decay of the critical field, 
for small values of this field. This type of behavior has been reported recently 
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in the experimental study of the copper doped transition metal dichalcogenide 
Cu x TiSe2 [14] and also in intercalated graphite compounds [15]. Since both 
systems potentially possess Dirac electrons as elementary excitations, it is con- 
ceivable that the peculiar properties of such electrons would be responsible for 
this common behavior of the critical field. In connection to this point, it would 
be extremely interesting to measure the critical magnetic field as a function 
of doping in Cu x TiSe2, in order to compare with the result derived from our 
B x Ar phase diagram. 



2 Model 

We investigate here the effect of applying a constant magnetic field Bz along 
the c-axis of a quasi-two-dimensional superconducting electronic system con- 
taining two Dirac points. Assuming that the active electrons correspond to 
these points, it follows that the electron kinematics will be described by a 
Dirac equation [1]. The electron creation operator will be V4,o-,a; where i — 1,2 
denotes the Dirac point and a —\, J,, the z-component of the spin, a = 1, N 
is an extra label, identifying the plane to which the electron belongs. Alterna- 
tively, in a multi-band system, a could be used to specify the electron band. 

The Lagrangian describing the system, in the presence of the external magnetic 
field is given by 



where A t is the vector potential corresponding to B, a are Pauli matrices and 
Hb is the Bohr magneton. The second and third terms, respectively, contain 
the coupling of the magnetic field to orbital and spin degrees of freedom. As in 
[3], we assume there is an effective superconducting interaction whose origin 
will not influence the results of this work, g is the superconducting coupling 
constant, which is supposed to depend on some external control parameter. 
In order to make the lagrangian smooth, we define g = X/N. We use the same 
convention for the Dirac matrices as in [3]. 

Performing a Hubbard- Stratonovich transformation we arrive at 




(1) 



C[*,a]=-- <T*<T + & a A* a 



(2) 
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where 



t d -d- 



A 



\ 



—d + d a 
a* d a, 
\ a* d- do, J 



(3) 



with do = i (^<9o + HbB), d± = i vf (hd± + i(e/c)A±) and <9-t = <9 2 ±i d\. The 
fermions are in the form of a Nambu field ^\ = (ip\^ a V4-|-a V'lj.a V4j.a) an d the 
auxiliary Hubbard-Stratonovitch field a satisfies the equation, 

<T=-g (>2|a ^lTa + ^l|a ^2Ta) • (4) 

This shows that is a Cooper pair creation operator. 

Integrating on the fermion fields, we obtain the effective action for a, namely 



S cS (\a\,B) = J d 3 x (-yM 2 ) -iAHnDet 



A [a, B] 



_A[a = 0,B = 0] 



(5) 



3 Effect of a magnetic field on the superconducting quantum phase 
transition at T = 



The effective potential per N corresponding to (5) may be obtained by a saddle 
point procedure, in which the field configurations assume their ground state 
average values that correspond to the classic lowest energy configurations. 

Observing that the determinant of the matrix A is given by 



det .4. = 



d - d+d- + \a\ 



(6) 



and choosing the asymmetric gauge, in which A = B(0,x), we infer from the 
effective action (5) that the effective potential will be 



e Vl '' ' A J (2tt) 2 J 2vr \ 



det A [a, B] 



det A [a = 0, B = 0] j ' 



(7) 
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where (6), in momentum-frequency space, is given by 



det A [a, B] = Ufwj + fi B B) 2 + h 2 v 2 F 



k 2 x +(k y + e -B{x) 



\ 2 

|2 



+ <t z + Bk 



(8) 

with k = Vph(e/c). 

In this expression, by a, we mean (0|<r|0). Analogously (x) is the average of 
the rc-coordinate in the lowest energy state of the relativistic Landau problem 
[16]. 

The effective potential in the presence of a magnetic field may be obtained by 
decomposing the logarithm in two parts and performing the shift of integra- 
tion variables k y + (e/c)B(x) — > k y and fkv + hbB — > Jtxj, in the first term. 
Introducing the momentum cutoff A/vf, we obtain, up to a constant, 



V M , B) JfL-^ ? + l {W? + BK f, ( 9 ) 
where a = 2nv F . 

The divergence can be eliminated as usual by the renormalization 



d 2 V cS 



dada* 



ii An A O. A r 



where is the renormalized coupling parameter, A c = 2a/3cr and 

f(B,a )= , ^ . (11) 
1 + ^ 

In the above equations, ctq is an arbitrary scale, the renormalization point. In 
the next section, we perform a renormalization group analysis, which shows 
that physical quantities do not depend on a . 

Using (10) in (9), we obtain the renormalized effective potential 

VMB) = K.fi^A wn l H + BK f (12) 



We may now determine the T = phase diagram as a function of the magnetic 
field by analyzing the minima of the effective potential above. For this purpose, 
we consider the derivatives of with respect to |<r|, namely 



(13) 
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and 

f( B ^o) , l rrr, , D . \ , / 2|<r 




(14) 



The minima of the effective potential will occur for the solutions of V^ s R (\a\,B) 
0, R (\a\, B) > 0. We, henceforth call these solutions A. Observe that 
A = |(0|<t|0)| and is, therefore, a superconducting order parameter. The pe- 
culiarities of the superconducting transition for T ^ in 2D are well-known 
and have been discussed extensively in the literature [17]. 

From (14) we see that a solution A = exists for 



X R < X C (B) = X c _ v \ +B _ , (15) 

where B = B(n/al). 

A solution with A ^ will only occur for Vl S R (A,B) = . From (13), we 
find 




r 2 pir^-^i ~ Bk - (16) 



Since the first term of (14) vanishes at this solution, we can readily see that 
V^b R (Ao, B) > 0, provided A is real, namely, for Xr > X C (B), where X C (B) 
is given by (15). We see that for this range of the renormalized coupling 
parameter, the superconducting gap (16) is a true minimum of the effective 
potential. 

The conclusion is that a quantum phase transition connecting a normal and 
a superconducting phase occurs at the magnetic field dependent quantum 
critical point X C (B), given by (15). Notice that in the limit 5^0 both the 
quantum critical point and the superconducting gap reduce to the ones found 
previously [3] in the absence of a magnetic field. Conversely, for each value of 
the physical coupling parameter Xr, there is a critical magnetic field above 
which superconductivity is destroyed. This observation allows us to infer the 
zero temperature phase diagram of the system, which is depicted in Fig. 1. 

Recent advanced techniques of controlled intercalation of Cu in TiSe 2 enabled 
the obtainment of the experimental measure of important physical parame- 
ters as a function of doping in this layered dichalcogenide [14]. It would be 
interesting to compare this phase diagram with corresponding experimental 
results of the critical field as a function of doping in Cu x TiSei. 
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We have shown in [3], that the mean- field phase structure obtained at zero 
magnetic field is robust against quantum fluctuations. The same arguments 
apply here and, therefore, we reach the same conclusion in the presence of an 
applied magnetic field. 



4 Renormalization Group Analysis 



In order to eliminate the high-momentum divergence, we have renormalized 
the theory by performing the subtraction (10) at the arbitrary scale <t . In 
this section, we show that the physical quantities obey renormalization group 
equations which show that they are actually independent of cr . 

By inserting the expressions of A c and of f(B,a ), in (12), we can show, 
for instance, that the effective potential satisfies the renormalization group 
equation 

( a "t a + ^ ir « = °- (17 » 

where the /9-function, defined by 

P = ^o^, (18) 



is given by 



where B = B(n/al). 



Ac (l + B) 2 



In analogous fashion, we can show that the superconducting gap (16) also 
satisfies the renormalization group equation (17), being therefore independent 

Of (Jo- 



Finally, by integrating the /5-function equation (18), we obtain the result 
1 f(a' ,B) _ 1 /«,£) 



\ R (a' ) \ c (a' ,B) A fl K) A C K,S)' 
for arbitrary scales a' and cr '. This explicitly shows that the combination 

1 f((Jo,B) 



(20) 



Ai? ((7 ) A c ((T , B) 



(21) 



does not depend on the renormalization scale a . From (16), then we explicitly 
see the scale independence of the superconducting gap A . 
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We conclude that the physical properties of the system will be determined by 
the physical coupling parameter Ar, which should be an experimental input. 
The value of the zero-magnetic-field quantum critical point A c must also be 
determined experimentally. Subsequently, the magnetic field dependence of 
the quantum critical point may be obtained from (15). 



5 Effect of a magnetic field on the superconducting phase transi- 
tion at T ^ 



Let us consider now the effect of an external magnetic field in the supercon- 
ducting properties of a quasi-two-dimensional system of Dirac electrons for 
T 7^ 0. In this case, the effective potential corresponding to (5) is given by 



I _ 1 2 „ j2l, oo 

V«(\<t\,B) = U--2 — -j J2 Un\(hw n + fi B B) 2 + h 2 v F \k\ 2 + \a\ 2 + BK 
A J (2tt) n= __ 



-In [hw 2 n + v 2 F \k\ 2 ]}, 



where u n are fermionic Matsubara frequencies, corresponding to the functional 
integration over the electron field. Observe that now we may no longer shift 
away the magnetic field coupled to the electron spin. This appears summed 
to the Matsubara frequencies. 

The phase diagram at T ^ is obtained from the minima of (22). For the 
purpose of determining these minima, we consider the necessary condition 
VIq (|<j|, B) = 0. Taking the derivative of (22) with respect to \a\ and perform- 
ing the Matsubara sum, we get 



v:„(w\,b) 



d 2 k 1 



(2tt) 2 E 



sinh ((3E) 



cosh (13 E) + cosh ((3 ji B B) 



0, 



(23) 



where E 



ilk 2 + \a\ 2 + Bk. 



Let us look for a superconducting phase. In such a phase, a nonzero solution 
of (23) (which we call A) is required. For this, the quantity between round 
brackets in (23) must vanish. After a change of variables, this condition leads 
to the equation for the superconducting gap A(T): 



1 = 



A /-a 



dE 



sinh (j3 E) 



a Jva 2 +Bk cosh (j3E) + cosh ((3 HbB) ' 
where A/vf is the high-momentum cutoff. 



(24) 



(22) 
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The second derivative of the effective potential evaluated at the solution of 
(24) is given by 



Kff (A, B) = 



2A 2 



a/VA 2 + kB 



sinh /VA 2 + kB 



cosh /VA 2 + kB + cosh (3n B B 



> 0. 



(25) 



This guarantees that the solution of (24) is indeed a minimum of the effective 
potential. 

In order to solve (24), we perform the integration in E. This may be done 
exactly and the result depends on the cutoff A. It is a well known fact that 
the inclusion of a finite temperature does not change the divergence structure 
of the theory. Indeed, the cutoff may be eliminated precisely by the same 
renormalization operation (10), which we used at T = 0. After this, we obtain 
the following expression for the superconducting gap A(T,B), as a function 
of the temperature and of the magnetic field, 



A 2 (T,B) 



k B T cosh -1 



cosh 



(HbB\ 
\k B Tj 



Bk, 



(26) 



where k B is the Boltzmann's constant and Ao = A(T = 0, B) is the zero 
temperature gap, given by (16). It is not difficult to show from the above 
equation that indeed, for T — > 0, we have A(T) — > A , as it should. 

Now, from the gap expression we may extract the critical temperature for the 
superconducting transition, T c . Indeed, using the fact that A(T C ) = in (26), 
we get, 



k B T c 



A 2 



Bk 



ln { 2 [ cosh (S)+ cosh (ff*)]} 



(27) 



Notice that for B = 0, we have A (0) / k B T c (0) = 2 In 2 in agreement with the 
result obtained in [3]. 



We may then re-express the gap as a function of T and T c as 



A 2 (T,B) = \ k B T cosh' 1 {2^- l) 



cosh 



'kB 



k B T c 



+ cosh 



\k B T c ) 



T c /T 



cosh 
(28) 



( »bB \ 
V k B T ) 



from which we confirm that the gap vanishes at T = T c 



In order to display the superconducting g £iP CIS db function of the temperature 
for different values of the magnetic field, we insert the expression for A , 
given by (16), in (26). Our results for Ar/A c = 2 are shown in Fig. 2 (the 
numerical results we present were obtained using v \ = 1.69 x 10 s (m/s) 2 and 
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a = (8/3) In 2 £; B T C (0) , which yields T c (0) = 2 K for \ R /\ C = 2 in the absence 
of magnetic field) . 



Moreover, by replacing (16) in (27) we may obtain T c self-consistently as a 
function of the coupling parameter. Our results are depicted in Fig. 3, which 
is the Tx Ar phase diagram of the system for different values of the magnetic 
field. As the magnetic field increases, the quantum critical point is shifted 
to the right. Our results suggest that this behavior may be experimentally 
verified for the copper-doped dichalcogenide Cu x TiSe2 [14]. 

From (27), we can also obtain the critical magnetic field as a function of the 
coupling parameter for different values of temperature, namely, the B x Ar 
phase diagram, which is shown in Fig. 4. We see that for very low temperatures, 
it reduces to the same phase diagram displayed in Fig. 1. However, as the 
temperature is increased, the superconducting region, which corresponds to 
the area below B c , becomes smaller. 

From (27), we may also obtain the B x T phase diagram for the quasi-two- 
dimensional superconducting Dirac electronic system. This is represented in 
Fig. 5. Particularly interesting is the linear behavior of the critical magnetic 
field for B > 0. We may derive explicitly from (27) the following expression 
for the critical field, in the small B region: 

* m „^^ (0) (!__*_), (29 ) 

where 




and 



A = 1 -4l^ 1 -^l- < 30) 




B C (T) exhibits a linear behavior of the critical field, which is indicated by the 
dotted lines for different values of the dimensionless parameter x = Xr/\ c in 
Fig. 5 (this is actually valid for \ R < 13 A c , when A is positive ). This differs 
from the quadratic behavior predicted by BCS theory. 

A linear decay of the critical field with the temperature similar to the one 
obtained here has been experimentally observed in intercalated graphite com- 
pounds [15] and also in the copper-doped dichalcogenide Cu x TiSe 2 [14]. Since 
graphene and also the transition metal dichalcogenides are well-known to pos- 
sess Dirac electrons in their spectrum of excitations, one is naturally led to 
wonder whether the presence of such electrons could explain such a behavior 
of the critical field. 
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It is interesting to observe that in expressions (26), (27) and (28) for the gap 
and T c , we can trace back the contributions from the spin and orbital couplings 
of the external magnetic field. These are given, respectively by the //# and k 
proportional terms. 



6 Conclusion 

Starting from the explicit expression for the superconducting gap as a function 
of temperature, magnetic field and of the effective superconducting coupling 
parameter, we have derived three phase diagrams for a superconducting quasi- 
two-dimensional system of Dirac electrons. 

The Tx Ar diagram displays a transition temperature with the same qualita- 
tive behavior as the critical temperature T c (as a function of doping) observed 
in high-T c cuprates in the underdoped region. It is conceivable that the effec- 
tive coupling parameter Ar, controlling the magnitude of the superconducting 
interaction, could be effectively determined by the amount of doping. In this 
case a direct comparison between our curve and the T c x doping experimen- 
tal results for the cuprates would be possible. Since the coupling magnitude 
should increase with doping, the curve of T c as a function of doping would 
have the same qualitative form as a function of the coupling Ar, wich by its 
turn would qualitatively agree with ours. The qualitative agreement with the 
cuprates data might be an indication that Dirac electrons play an important 
role in the mechanism of high-T c superconductivity. It would be interesting 
to investigate the magnetic field dependence of the superconducting dome in 
cuprates and especially of the T = critical doping determining the onset of 
superconductivity. 

The B x T phase diagram, by its turn shows a linear decay of the critical field 
as a function of temperature, for small fields, which differs from the corre- 
sponding quadractic behavior predicted by BCS theory. It is quite interesting 
that two materials that are supposed to have Dirac electrons as elementary 
excitations, namely, the copper-doped dichalcogenide Cu x TiSt2 [14] and in- 
tercalated graphite [15], both present the same type of linear behavior of the 
upper critical field as a function of temperature, in the small field regime. This 
might indicate that Dirac electrons are the common cause of this behavior in 
both materials. We intend to explore this point more profoundly in a next 
publication. 

Finally, the B x Ar phase diagram presents a quantum critical line, corre- 
sponding to a magnetic field dependent quantum critical parameter X C (B), 
connecting the normal and superconducting phases. Again, we may assume a 
direct connection between the coupling parameter Ar and the doping parame- 
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ter, for instance, in the case of Cu x TiSe2- In this case, we would actually have 
in Fig. 4 a phase diagram B x x, where x is the doping parameter. It would 
be extremely interesting to have experimental curves of the critical magnetic 
field as a function of doping in dichalcogenides such as Cu x TiSe2, in order to 
compare with our theoretical results. 
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Fig. 1. The zero temperature phase diagram of the system. 
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Fig. 2. The superconducting gap A/ks as a function of the temperature for several 
values of the magnetic field. 
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Fig. 3. T x Ar phase diagram. The superconducting critical temperature T c as a func- 
tion of the dimensionless coupling parameter Ar/A c for several values of magnetic 
field. 
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Fig. 4. B x Ar phase diagram. The critical magnetic field B c as a function of the 
dimensionless coupling parameter Ar/A c for several values of the temperature. 
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Fig. 5. B x T phase diagram. The critical magnetic field B c as a function of the 
temperature for several values of the dimensionless coupling parameter x = Ar/A c . 
The dotted lines in the figure indicate the linear behavior of B given by (29) as 
B -> 0. 
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